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ABSTRACT 

We present both photometry and spectra of the individual images of the quadruple 
gravitational lens system Q2237+0305. Comparison of spectra obtained at two epochs, 
separated by ~ 3 years, shows evidence for significant changes in the emission line to 
continuum ratio of the strong ultraviolet CIV A1549, CHI] A1909 and Mgll A2798 lines. 
The short, ~ 1 day, light -travel time differences between the sight lines to the four 
individual quasar images rule out any explanation based on intrinsic variability of the 
source. The spectroscopic differences thus represent direct detection of microlensing- 
induced spectroscopic differences in a quasar. The observations allow constraints to be 
placed on the relative spatial scales in the nucleus of the quasar, with the ultra-violet 
continuum arising in a region of Sa 0.05 pc in extent, while the broad emission line 
material is distributed on scales much greater than this. 
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1 INTRODUCTION 

Fluctuations in the brightness of cosmologically distant 
sources due to the action of gravitational microl ensing is 
now a well established observ ational phenomenon (Irwin et 
|al. 1989 ; Corrigan et al. 1991), also occurring, spectacularly, 



2 MICROLENSING AND Q2237+0305 



Discovered during the CfA redshift survey (Huchra et al 



1985) the Q2237+0305 system consists of four images of a 
back ground quasar, z ~ 1.69, separated by ^ 2".0 (Yec 



within our own Galaxy (cf Alcock et al 1996). Utilizing ob- 
servations of the quadruple gravitational lens Q2237+0305 
we present evidence of another manifestation of microlens- 
ing; modification of the observed spectral characteristics of 
a high-redshift quasar due to differential amplification of 
source structures with different spatial scales by objects of 
stellar mass along the line-of-sight. 

Section [| presents an overview of microlensing and 
the properties of the quadruple lens Q2237+0305 that 
make it an almost ideal target for microlensing investiga- 
tions. The spectroscopic and photometric observations of 
the Q2237+0305 system are presented in Section 3, with a 
discussion of the photometry appearing Section B The re- 
sults of the spectral extraction are presented in Section 
and an interpretation follows in Section ^| The basic data re- 
duction and the relatively complex procedure for extracting 
the spectra of the individual quasar components is presented 
in an appendix to this paper. 
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The images are centred on the core of the lensing 
galaxy, a bright nearby, z ~ 0.04, barred spiral. The galaxy 
light in the inner regions of the Q2237+0305 syste m is dom i- 
nated by stars in the central bulge of the galaxy ( Yee 1988 ), 
and the galaxy surface brightness profile was modeled us- 
ing a de Vaucoulers (1953) r* profile, wi th a core rad ius of 
r e = 7"0, and an ellipticity of e = 0.31 (iRacine 199lh. The 



stellar density along the sight-lines to the four quasar images 
is considerable and microlensing by stars within the galaxy is 
expected to significantly perturb the observed image fluxes. 

The microlensing scale-length for an isolated star of 
mass M, its Einstein radius rjo, is defined to be; 



4GM D os Di, 



D 



(1) 



where Dij are angular diameter distances between the ob- 



server (o), lens (1) and source (s) (Schneider et al. 1992). 
This length is defined in the "source plane" at a distance 
D os from an observer. A point-like source passing within 
this projected radius of a microlensing mass is amplified by 
a factor of at least 1.34, although the angular scale of mul- 
tiple imaging induced by such lenses is ~ 10~ 6 ", far below 
present detection limits with optical telescopes. In practice, 
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therefore, observations of microlensed systems are confined 
to the monitoring the image brightness fluctuations. 

The light curve induced by the passage of an isolated 
microlens possesses a very simple form, with a characteristic 
time- scale, r: 



V cff : 



(2) 



where V e ff is the effective veloc ity of the source across the 
source plane (Kayser et al. 1986). In such simple cas es r can 



be use d to c haracterize the mass of the microlenses (Wamb 
sganss 992). 

In a high optical depth regime, as is the case for 
Q2237+0305, the action of the individual lensing masses 
combine in a highly non-linear fashion, and the time-scale of 
individual "events" no longer reflects the time-scale given by 
Equation ^. Instead, the light curve of a background source 
exhibits complex variability behaviour, including asymmet 



ric fluc tuat ions (Paczyriski 1986; Wambsganss 1990 
et al. lQ 



Lewis 



Notwithstanding the complexities arising from the 
high optical depth of stars through the lensing galaxy, 
Q2237+0305 is unique among known gravitational lens sys- 
tems in that the lens is located very close to the observer 
and thus the light-travel time differences to the four im- 
ages are only ~ 1 day. Intrinsic variations in the luminosity 
of the quasar thus manifest themselves in all four images 
within ~ 1 day and differential variability between images 
with time-scales ^> lday reflects the effects of microlensing. 
It is this effective decoupling of the observed time-scales for 
intrinsic variability and microlensing by stars in the galaxy 
that allowed photometric monitoring of Q2237+ 0305 to pro 



duc e th e first detection of a microlensing event ( [rwin et al 
198fi|)~l 

Chang (1984) demonstrated that the amplification of 
a source is dependent on its size relative to the Einstein 
radius of the characteristic lensing mass. At caustic crossings 
a point source is amplified by an infinite amount, but any 
physical extent leads to a finite amplification. Sources whose 
scale size is only a fraction of an Einstein radius can be 
amplified by large factors, while sources whose scale is much 
greater than an Einstein radius are amplified by a negligible 
amount. 

For Q2237+0305 the Einstein radius (Equation [j]) for a 
star of Solar mass, projected into the source plane, is 



7/0 ~ 0.05/i 5O ' pc. 



(3) 



Through-out we assume a standard Freidmann-Walker Uni- 
verse, with an Q = 1 of matter distributed smoothly. The 
cosmological constant, A, is assumed to be zero. Within the 
framewo rk of the "s tandard model" for the central regions of 
quasars (Rees 1984) this microlensing scale is smaller than 
the extent of the broad emission line region, 0.1-1.0 pc, 
but significantly in excess of the size of the ultraviolet- 
optical continuum-producing region, < 10" 3 pc. Thus, dif- 
ferential microlensing amplification of the quasar continuum 
and emission line producing regions by objects of stellar 
mass within the lensing galaxy of the Q2237+0305 system 
is expected. This should be observable as a time-dependent 
variation in the relative strengths of the broad emission 
lines and the continuum, and should also be coupled with 
microlensing-induced photometric variability. Such observa- 



tions would provide verification of the microlensing hypoth- 
esis and provide a powerful diagnostic of quasar structure 
on extremely small scales. 



3 OBSERVATIONS 

Quasi-simultaneous spectroscopic and imaging observations 
of Q2237+0305 were obtained at two epochs, separated by 
approximately three years, using the 4.2m William Herschel 
Telescope (WHT) at the Roque de los Muchachos Observa- 
tory, La Palma. 



3.1 CCD Imaging 

Direct images of Q2237+0305 were obtained at both epochs 
through a KPNO broad-band R filter using the auxiliary 
port at the Cassegrain focus of the WHT. On 1991 August 
12, between 04:00 and 04:40 UT, six 300 s exposures us- 
ing an EEV CCD were acquired. The seeing was excellent 
with a measured FWHM = 0"5. The second epoch images 
consisted of a 300 s exposure obtained on 1994 August 14, 
between 23:50 and 23:54 UT, also employing an EEV CCD. 
Again, atmospheric conditions were excellent with the seeing 
measured to have a FWHM = 0"4. The image scale with the 
EEV CCD at the auxiliary port of the WHT is 0" 11 /pixel. 
Bias frames, and an exposure of a region of sky devoid of 
bright stars for flat-fielding were obtained during evening 
twilights. Bias subtraction and flat-fielding procedures were 
undertaken using IRAF routines. 



3.2 Spectroscopy: Instrumental Configuration 

The spectroscopic observations were acquired w ith the 
ISIS double-beam spectrograph (Clegg et al. 1992) at the 
Cassegrain focus of the WHT. On 1991 August 14 the spec- 
trograph was configured with the 5400 dichroic, splitting the 
red and blue beams at ~ 5400A. The R158B grating in the 
blue arm combined with a thick 1152 x 1242 pixel EEV CCD 
with 22.5/jm pixels produced a dispersion of 2.70 A/pixel. 
The grating angle was set to give a central wavelength of 
~ 4500 A. The wavelength coverage in the blue is truncated 
at ~ 3500A by the combined effects of atmospheric absorp- 
tion and the rapidly falling sensitivity of the CCD, and in 
the red, at ~ 5500A, by the dichroic cross-over wavelength. 
The R316R grating in the red arm combined with a second 
thick 1152 x 1242 pixel EEV CCD with 22.5/^m pixels pro- 
duced a dispersion of 1.40 A/pixel. The central wavelength 
was set at ~ 7800A, giving a wavelength coverage AA7000- 
8600 A. The spatial scale on the detector along the slit was 
0"335/pixel for both red and blue arms. 

The 1994 August 17 observations were obtained with a 
very similar instrumental configuration. Both red and blue 
gratings, central wavelengths and wavelength coverage were 
as for the 1991 observations. The 5700 dichroic split the red 
and blue beams at ~ 5700A, changing the wavelength cov- 
erage between the two arms slightly. The blue-arm detector 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by AURA under co-operative agree- 
ment with the National Science Foundation. 
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Figure 1. A grey-scale representation of the R-band WHT images of Q2237+0305 from the 1991 and 1994 epochs. The images have 
identical scales and orientations, and show a region ~ 3"25 on a side with north to the top and east to the left. The legend in the 
lower left-hand corner of the August 1991 images shows the image designations of Yee (1988). Image A can be seen to have brightened 
substantially over the three year period. 



had been upgraded to a 1024 x 1024 pixel thinned TEK CCD 
with 24fim pixels, decreasing the dispersion to 2.88 A/pixel 
but producing much improved sensitivity. The spatial scale 
on the detectors along the slit was 0"335/pixel in the red 
arm and 0"358 in the blue arm. 

Observations at both epochs therefore covered wave- 
lengths AA3500 - 5400 and AA7000 - 8600, which, at the 
redshift of z = 1.69 for Q2237+0305, include the prominent 
rest-frame emission lines of CIV A1549, CIII] A1909 and 
Mgll A2798 at wavelengths of ~ 4175A, 5145A and 7540A, 
respectively. Unfortunately, the strong telluric absorption 
band at ~ 7600A coincides with the red wing of the Mgll 
emission line. 



3.3 Spectroscopy: Observations 

The atmospheric conditions during the spectroscopic obser- 
vations at both epochs (1991 August 15 00:00-05:00UT and 
1994 August 17 05:00-05:30UT), were excellent. The seeing 
derived from the spatial extent of the component quasar 
images was measured to be < 0"8 at both epochs. A slit- 
width of 0"70 was used in 1991, producing spectral resolu- 
tions (FWHM of an unresolved arc line) of 8.1 A and 4.1 A in 
the blue and the red arms respectively. In 1994 a slit- width of 
0"91 was employed giving spectral resolutions of 5.7A, and 
4.lA in the blue and the red arms respectively. Q2237+0305 
was observed either side of the meridian with spectra ob- 
tained at airmasses of 1.1 to 1.4. The lower resolution in the 
blue-arm during the 1991 observations is probably a result 
of non-optimal spectrograph focus. 

Spectra of pairs of images on the same side of the galaxy 
nucleus were obtained in order to minimize contamination 
from the galaxy core. Observations of sufficient pairs were 
acquired to ensure some redundancy in the observations of 
each component. The longer 1991 series comprised image 
pairs A+D, B+C, B+D, and A+C while the short sequence 
in 1994 (obtained as a Service allocation) comprised A+C 



Year 


Image Pair 


Slit PA 


Exposure 


Airmass 


Seeing 


1991 


A+D 


302.8 


1800 


1.386 


0'(80 


1991 


A+D 


302.8 


1800 


1.254 


0'{80 


1991 


B+C 


290.8 


1800 


1.144 


0':80 


1991 


B+C 


290.8 


1800 


1.125 


0'!80 


1991 


B+D 


8.8 


1800 


1.111 


0':80 


1991 


B+D 


8.8 


1800 


1.134 


0".80 


1991 


A+C 


27.8 


1800 


1.218 


o':so 


1991 


A+C 


27.8 


1800 


1.316 


0".80 


1994 


A+C 


27.8 


1500 


1.109 


0'(65 


1994 


B+C 


290.8 


1500 


1.108 


0".65 



Table 1. Journal of spectroscopic observations of Q2237+0305 
taken on the nights of 1991 August 14/15 and 1994 August 17/18. 



and B+C. The four individual quasar images were clearly 
visible on the acquisition TV and the spectrograph slit, ori- 
ented to the appropriate position angle for the particular- 
image pair, was aligned across each image pair using visual 
inspection. A log of the spectroscopic observations is given in 
Table [l} It was necessarily impossible to acquire the spec- 
tra of the image pairs with the spectrograph slit aligned 
along the parallactic angle. Differential atmospheric disper- 
sion was thus expected to produce a systematic loss of light 
from the slit at the bluest wavelengths. 

Bias frames and exposures of quartz and tungsten lamps 
were o btained at the start of each ni ght and flux stan- 
dards (Massey et al. 1988; Oke 1990) were observed 



at 



the end of the night. Exposures of Cu+Ar (blue arm) and 
Cu+Ar/Cu+Ne (red arm) lamps were taken following each 
exposure of an image pair or a standard star. 
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Seeing 


A 


B 


C 


D 


1991 
1994 


0"A 


17.53 
17.12 


17.23 
17.20 


18.11 
18.00 


18.28 
18.23 



Table 2. The R-band image magnitudes of Q2237+0305 at the 
1991 and 1994 epochs. The errors in these measurements is 0.05 
mag. 



4 PHOTOMETRY 

Figure |l] shows a grey-scale representation of the reduced 
R-band frames from the 1991 and 1994 epochs. Relative 
photometry of the individual quasar images was performed 
using the procedures described in Corrigan (1993). Expo- 
sures of standard stars were not obtained at either epoch 
and the frames were calibrated using stars visible both in 
frames from Corrigan (1993) and the new images. The un- 
certainty in the calibration is estimated to be ±0.05 mag. 
R-band magnitudes for the four quasar images on the R- 
band system used in Corrigan et al (1991), are given in Ta- 
ble ^. The significant difference in the relative brightness of 
images A and B between the two epochs is evident from in- 
spection of Figure [l], and the change in relative brightness, 
ArriA-B = 0.38 ± 0.07 mag is well determined, independent 
of the calibration procedure. 

Relative, and calibrated, magnitudes for the images at 
both epochs are consistent with the photometric time se- 
ries (27 epochs) for the period 1990 to 1993, obtained by 
0stensen (1994) at the Nordic Optical Telescope, which 
shows image A brightening to become comparable to im- 
age B at the end of 1992. 



5 SPECTROSCOPIC REDUCTIONS 

The relatively complex nature of the data, consisting of two 
merging quasar spectra superimposed upon an extended 
galaxy background, precluded the application of standard 
data reduction software to extract the spectra of the in- 
dividual quasar components. A new reduction technique, 
therefore, was developed for application to this data. This 
technique is presented, in some detail, in Appendix 



5.1 Spectral Fitting 

The image separations, seeing and (x 2 )i the average 
"goodness-of-fit" value (Appendix A4), from the model- 



fitting procedure applied to the blue and red spectroscopic 
observations of each image pair are presented in Table ^| 
The seeing is defined as the FWHM of the Gaussian core 
of t he pro file model fitted to the two quasar spectra (Sec- 
tion A3.1). For comparison, image separations derived from 



Hubble Space Telescope imaging of the system (Rix et al 



1992) are presented in Table |3| 

We note that the red spectra in Table |^ possess system- 
atically greater estimates of the width of the seeing profile 
when compared to the equivalent blue spectra. This may be 
due to an under-estimate of the galaxy component in the 
wings of the quasar profile, driving profile width to larger 
values. This, however, was not readily apparent in "by-eye" 



inspections of the fits as a function of wavelength. Such 
model dependencies will not effect the over-all conclusions 
presented in this paper. 



5.2 The Spectra 

As an illustrative example of the extracted spectra, Figure ^ 
presents all four blue spectra of image A obtained in 1991. 
The spectra are normalized such that the average flux in the 
continuum region 4400 — > 4900A is the same. To within the 
noise, there is excellent agreement between the independent 
observations of the spectrum, with typical differences being 
~ 10%. 

It should be noted that the spectra show no apparent 
distortion due to the effects of atmospheric dispersion. The 
A+D and A+C image pairs were observed at very similar 
zenith distances, but either side of the upper culmination. 
At these air masses the resulting slit orientations lay within 
~ 15° of the parallactic angle, resulting similar atmospheric 
dispersion for each spectrum. This was not the case for the 
1994 observations, both of which were observed very near the 
upper culmination. Atmospheric dispersion, coupled with a 
narrow slit and good seeing, produced significant loss of light 
at the blue end of the spectra. 



5.3 Scaled Spectra 

The narrow slits and the requirement to orient the position 
angle away from the parallactic, (Section^, mean that spec- 
tra of the individual images cannot be compared using an 
absolute flux scale. However, the principal effect predicted 
to arise from differential microlensing of the continuum and 
emission line producing regions will manifest itself as differ- 
ences in emission line equivalent -width among the spectra 
of the images. 

Any such differences in equivalent-widths between the 
spectra should be evident visually if the spectra are scaled 
such that their continuum levels are coincident. Equivalent- 
width differences will then manifest themselves as emission 
line strength changes. Figures ^|-^| show the 1991 epoch 
blue and red spectra of such continuum-scaled spectra for 
the image pairs, A+C, A+D, B+C and B+D respectively. 
The spectra are plotted following the multiplicative scal- 
ing of the fainter spectrum such that the mean fluxes in 
the "continuum" regions 4650 — ► 4900 A (in the blue), and 
7000 — > 7080A in the red spectra are equal. 

Spectra of the image pair A+B (Figure ^) show no obvi- 
ous systematic difference in emission line strength. However, 
Figure ^ shows image D has significantly stronger emission 
lines than image A and similarly Figure H demonstrates that 
image C has stronger emission lines than image B. Consis- 
tency in the emission line strength differences between the 
images is demonstrated in Figure where image D is seen to 
possess a very marked excess of emission compared to the 
spectrum of image B. In addition to the three prominent 
emission lines of CIV, CM] and Mgll, evidence for an excess 
at 4455 A, coincident with Hell A 1640 is present. 
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Blue Spectra Red Spectra 

Images Year A FWHM (x 5 ) A FWHM (x 5 ) A H st 



A+D 


1991 


1".07 ± 0^25 


1".06 ± 0".U 


1.23 


1^07 + 0^23 


1^13 ± 0^14 


1.91 


1"01 


A+D 


1991 


1^09 ± 0".26 


0'(94±0'(11 


1.29 


1"05 ± 0'!24 


1"11 ±0'(13 


1.91 


1"01 


B+C 


1991 


1".38 ± 0".22 


0".83 ± 0^13 


1.34 


1"44±0"19 


1'iOl ± O'ill 


1.48 


1"39 


B+C 


1991 


1".36 ± 0".25 


0'(90 ± 0"18 


1.40 


1"A2 ± 0'!22 


l'(05 ± 0^13 


1.74 


l".39 


B+D 


1991 


1".23 ± 0".27 


0".72 ± 0'(05 


1.21 


1".21 ± 0".30 


l'(05 ± 0'(10 


1.77 




B+D 


1991 


1".20 ± 0".31 


0".72 ± 0".06 


1.25 










A+C 


1991 


l".33 ± 0".22 


0'(84 ± 0'(14 


1.40 


1".34 ± 0".20 


l'{12 ±0'(13 


1.60 


1"35 


A+C 


1991 


1".28 ± 0".22 


0"89 ± 0".ll 


1.30 


1".35 ± 0".21 


1"17±0'(13 


1.61 


1"35 


A+C 


1994 


1".33 ± 0".09 


0"57 ± 0'(03 


1.46 


1".24 ± 0'!l2 


0"99 ± 0'(05 


1.60 


1"35 


B+C 


1994 


1".33 ± 0'(30 


0".79 ± 0".W 


1.06 


l'{34 ± 0'(13 


0'(90 ± 0'(08 


1.47 


l'{39 



Table 3. Spatial quantities (image separation and seeing) and average goodness— of— fit derived from the model— fitting to each of the 
blue and red spectroscopic observations. The one-sigma errors are derived from the total spectral fit. The final column presents the 
image separations, Ahsti from the HST observations of Rix et al. (1992). The quoted error in the image positions as derived from HST 
imaging is ^ 0"015. 




4000 4200 4400 4600 4800 5000 



Wavelength A 



Figure 2. The four blue spectra of Image A obtained in 1991. The spectra have been corrected for the instrument response and normalized 
such that the mean flux in the continuum region 4400 -> 4900A is identical. 



5.4 Emission Line Equivalent Widths 

The emission line equivalent width (EW) provides a quan- 
titative measure of emission line strength: 

EW= f^mx)-nx ) ] dX (4) 



where S (A) is the flux at wavelength Aand / (A) is the flux in 
the underlying continuum at a wavelength A. Wavelengths, 
Ai and A2 are chosen to bound the line of interest. 

Establishing the true continuum level in rest-frame ul- 
traviolet spectra of quasars is difficult due to the presence 
of weak emission features [cf the composite quasar spectrum 
of Francis et al. (1991)]. However, establishing the absolute 
equivalent width is not required, rather, a consistent mea- 
sure of equivalent width that allows the relative emission line 
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a 
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<0 
(h 

X 



CO 



d 




d 



V 

<D o 
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d 
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Blue Arm 



_i_ 



_i_ 



Red Arm 




_i_ 



_i_ 



4500 



5000 7400 
Wavelength (A) 



7600 



7800 



Figure 3. The blue (left panel) and red (right panel) ISIS spectra for images A + C from the 1991 observations. In the top panels the 
spectrum of C has been scaled to match the spectrum of image A in the (4650 — > 4900A) in the blue spectrum, and ^7000 — > 7080A) 
in the red spectrum. The brighter A image is presented as a thick line, where as the scaled C spectrum is a thin line. The lower panel 
presents the difference in these spectra, defined as the spectra of A minus the scaled spectrum of C. 



Line 


Xl,l 


Xl,2 




X u ,2 


CIV 


3970 


4040 


4290 


4300 


CIII] 


5020 


5060 


5220 


5260 


Mgll 


7450 


7480 


7700 


7739 



Table 4. For the emission lines, xi and x u represent the bounds 
on the regions used in the estimation of the continuum. The equiv- 
alent width is measured between x; 2 and x u i> except for Mgll 
where the summation was truncated at 7590A to exclude the at- 
mospheric A-band absorption. 



strength in the component spectra is all that is necessary. A 
continuum was estimated for each emission line based on the 
median flux within "continuum" bands on either side of the 
line. The continuum was set by the straight line connecting 
points defined by the median flux in each band at the mid- 
point of the wavelength range of each band. The EW was 
then calculated numerically according to Equation ^ with 
the summation extending from the long wavelength edge of 
the blue continuum band to the short wavelength edge of the 
red continuum band (Table ^). In the case of the Mgll line 
the summation was truncated at 7590A to avoid the strong 
telluric absorption feature at ~ 7600A. 

Table I presents equivalent width measurements of CVI, 
CIII] and Mgll for spectra taken in 1991 and 1994. The errors 



Image 


Year 


CIV 




CIII] 






Mgll 




A 


1991 


126.11 ± 1.52 


43 


73 ±1 


34 


28 


65 ±0 


96 


D 


1991 


190.53 ± 4.35 


61 


28 ± 2 


96 


19 


83 ± 1 


92 


A 


1991 


132.23 ± 1.72 


38 


78 ±1 


12 


33 


10 ±0 


92 


D 


1991 


203.47 ±5.08 


64 


02 ± 2 


61 


13 


09 ± 1 


72 


B 


1991 


68.32 ±0.65 


20 


91 ± 


89 


18 


08 ± 


53 


C 


1991 


116.87 ±2.06 


47 


23 ± 2 


11 


28 


53 ± 1 


18 


B 


1991 


77.00 ±0.61 


25 


66 ± 


86 


19 


65 ± 


55 


C 


1991 


120.34 ± 1.93 


45 


66 ± 2 


30 


32 


61 ± 1 


17 


B 


1991 


69.59 ±0.68 


28 


96 ± 


93 


19 


26 ± 


64 


D 


1991 


148.88 ± 2.65 


83 


12 ±2 


11 


19 


23 ± 1 


92 


B 


1991 


69.70 ±0.84 


29 


03 ± 


91 








D 


1991 


179.62 ± 3.73 


93 


29 ± 3 


31 








A 


1991 


99.91 ± 1.03 


40 


61 ± 1 


06 


28 


05 ± 


73 


C 


1991 


112.41 ± 2.12 


36 


92 ± 1 


76 


25 


05 ± 1 


16 


A 


1991 


118.33 ± 1.33 


38 


20 ± 1 


16 


27 


48 ± 


77 


C 


1991 


83.37 ±2.20 


37 


45 ± 2 


01 


24 


76 ±1 


29 


A 


1994 


102.21 ± 0.43 


36 


44 ±0 


63 


25 


93 ± 


46 


C 


1994 


151.80 ± 1.60 


51 


73 ±1 


61 


29 


96 ± 


91 


B 


1994 


101.45 ± 1.12 


30 


51 ± 1 


64 


17 


36 ± 


58 


C 


1994 


141.34 ±3.18 


15 


51 ± 1 


93 


20 


61 ± 1 


10 



Table 5. Equivalent width measurements for images at both 
epochs. The error estimate is calculated assuming a continuum 
placement uncertainty of 2% for CIV and CIII, and a 5% error 
for Mgll. 
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Figure 4. As for Figure presenting the spectra for image A and image D in 1991, which has been scaled so that the continua match. 
The thinner line represents the scaled spectrum of the D image, while the thicker line is for image A. The difference between these 
spectra, presented in the lower panels, is defined to be the spectrum of the A image subtracted from the scaled D spectrum. 



Image 


Year 


CIV 


cm] 


Mgll 


A 


1991 


119.14 ± 12.15 


41.33 ± 3.80 


29.32 ± 2.22 


B 


1991 


71.15 ±3.42 


26.14 ±3.31 


19.00 ±0.67 


C 


1991 


108.33 ± 14.50 


41.82 ±4.67 


27.74 ±3.18 


D 


1991 


180.63 ±20.18 


75.43 ± 13.31 


47.38 ± 3.08 


A 


1994 


102.21 


36.44 


25.93 


B 


1994 


101.45 


30.51 


17.36 


C 


1994 


146.57 ±5.23 


48.62 ± 3.11 


25.28 ±4.67 



Table 6. Mean values of the emission line equivalent widths from 
independent spectra. The la errors represent the measured error 
in the mean. Only one spectrum of each of images A and B was 
obtained in 1994, and the errors in the values of the equivalent 
widths are assumed to be similar for those found for Image C. 



are calculated assuming an uncertainty in the continuum 
placement of 2% for CIV and CIII, and 5% for Mgll. Mean 
equivalent widths for each image at each epoch are given in 
Table |(| together with the la errors in the mean. 



5.5 B mioaion Line Ccntroid Shifta 



Nemiroff (1988) and Schneider and Wambsganss (1990) pre- 
dicted that gravitational microlensing of substructure in the 
broad line emitting region would result in significant profile 
differences between the same emission line in different im- 



ages. Selective microlensing enhancement of portions of the 
emission line region undergoing ordered rotation, or outflow, 
could result in shifts of ~ 1000km s _1 in the central wave- 
lengths of emission lines. 

Any differences in emission line centroids can be quanti- 
fied by cross-correlation. The cross-correlation between two 
discrete signals [ii and Vj is defined to be: 



N-l 

C(n,l>)k = ^2 Vk+jVj- 
3=0 



(5) 



Computation of the cross-correlation function using ob- 
served galaxy spectra and zero-redshift t emplates is a stan- 



dard technique in establishing redshifts (Tonry and Davie; 



197£). Here, emission lines in different images arising from 
the same transition were cross-correlated. The quasar con- 
tinuum was first subtracted from each emission line using 
the continuum definition procedure employed to perform the 
measurement of the equivalent widths. Then, the pairs of 
continuum subtracted emission lines were cross-correlated 



using the Numerical Recipes CORREL routine (Press et al 



198£). Due to the presence of the telluric A-band absorp- 
tion, only the region blue-ward of 7590A was used in the 
analysis of the Mgll velocity characteristics. 

A Gaussian function was fitted to the central portion of 
the resulting cross-correlation function in order to provide 
a measure of the peak location. The function took the form: 
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C(x) = Ae 



(6) 



where x c is the central position of the peak of the distri- 
bution of width c. A is a normalization factor, and d is a 
free parameter. The fitting algorithm was an implementa- 
tion of the Levenberg-Marquardt minimization method, as 
presented by Press et al. (1988) in MRQMIN. The peak-fitting 
was undertaken using apertures ranging from ten to eighty 
pixels either side of the peak. The location of the peak cen- 
tre depended somewhat upon the aperture used, as the form 
of the cross-correlation peak is not perfectly represented 
by Equation ^| ie large apertures tended to fit the cross- 
correlation peak wings, while small apertures didn't possess 
the curvature sufficient to accurately identify the peak. For 
these reasons various apertures were applied to each cross- 
correlation peak. This resulted in the selection of a separate 
aperture to to reproduce the form of each peak. 

The sensitivity of the cross-correlation analysis with 
spectra of the signal-to-noise ratio available is approxi- 
mately 50 km s _1 (la). The velocity shifts for any of the 
three emission lines between any of the image pairs were 
all consistent with no velocity offset. For each line in the 
1991 data, the mean velocity was; 4. 6± 20.9km s for CIV, 
48.1 ±88.9 km s _1 for CIII and 18.0 ±27.3 km s" 1 for Mgll. 
The uncertainties quoted are la and are based on the un- 
certainty of the peak-fitting procedure described above. No 
systematic offset was seen for any of the cross-correlation 
peaks and these values suggest that there are no significant 
differences in the central velocity of the lines. 

Consideration was also given to whether emission line 
profile differences were present between images. Such differ- 
ences are expected due to some preferential ampli fication of 



have a radio brightnes s which is compa rable to the other 
images in this system (Falco et al. 1996). The photometry 



any substructure in the broad emission l ine region ( Nemirofl 



1985 



Schneider and Wambsganss 199C). However, the lim- 



ited signal-to-noise ratio of the data precludes the derivation 
of astrophysically useful constraints on possible line shape 
differences. 



6 DISCUSSION 



6.1 Quasar Emission Regions 



Table ^ and the spectra shown in Section 5.3 demonstrate 
that during the 1991 observations, strong equivalent width 
differences existed between the spectra of the four images. 
This result is consistent with standard models of quasar 
structure, with the scale of the broad line emitting region 
significantly greater than that of the continuum source. 
Schemes, such as the star-burst model, where emission lines 
originate in narrow supernov ae shells embedded i n a large 
continuum emitting source ( |Cid- Fernandez 1995| ), cannot 
account for the systematic lack of enhancement of the line 
flux during microlensing. 

Table u\ presents the equivalent widths from both epochs 
of observation, normalized with respect to the equivalent of 
each line in the spectrum of image C. During the 1991 obser- 
vations image A had equivalent widths of order unity when 
compared to those of image C, while each value in image B 
was ~ 0.65 that of C. Image D, on the other had equivalent 
width measures of ~ 1.7 those of C, reflecting the microlens- 
ing deamplification of its continuum. This interpretation is 
consistent with recent observations which show image D to 



of system at the two epochs illustrate that image A had 
brightened over the three years, becoming as bright as im- 
age B. This is reflected in the values of Table ^. Essentially, 
the equivalent widths of image B, relative to those of im- 
age C, are unchanged between 1991 and 1994, while image 
A, on the other hand, has shown a strong decrease in the 
line-to-continuum strength in both CIV and CIII]. The EW 
measurements can be explained if image A underwent a mi- 
crolensing event between the two epochs. The measure of 
the Mgll equivalent width in image A, relative to that of 
image C, does not change between the 1991 and 1994 ob- 
servation. The reason for this constant value in Mgll is still 
unknown. 

The ratio of the equivalent widths presented in Table ^ 
indicate that, at the time of observation, the spectral slope of 
the quasar images was independent of the level of microlens- 
ing amplification, with an equal level of enhancement in the 
red and blue regions of the spectrum. In the rest frame of the 
quasar, the spectra presented here cover the range 1500A to 
2800A, which, if one adopts a temperature structure of the 
form T(R) oc as appropriate for the inner, radiation- 

dominated regi on of an optically thick, ge ometrically thin 
accretion disk (Shakura and Sunyaev 1973), implies that 

-R2800 ^ g 

-R1500 

For the corresponding effective temperatures, these emis- 
sion regions are £ 10 15 cm in extent, both the red and blue 
light comes from a fraction of an Einstein radius (Equa- 
tion [j]), and light curves in both colours should exhibit sim- 
ilar properties. The difference in scale size will only be ap- 
parent as colour differences when a caus tic passes over the 
source (Wambsganss and Paczynski 1991) and it is probable 
that all the sources lie in inter-caustic regions. 



6.2 Structure in the Broad Line Region 

To the limit of the sensitivity, ~ 150 km s _1 , no significant 
velocity differences between the same emission lines in the 
spectra of different images in the Q2237±0305 system were 
found. Thus, any substructure in the broad emission line 
region must be large enough to be insensitive to the effects 
of microlensing. With Equation ^, this implies substructure 

in the broad line region must be on scales J> 0.05/i BO 5 pc. 
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Image 


Year 


CIV 


CIII] 


Mgll 


Year 


CIV 


CIII] Mgll 


A 


1991 


1.10 


0.99 


1.06 


1994 


0.70 


0.75 1.02 


B 


1991 


0.66 


0.63 


0.68 


1994 


0.69 


0.63 0.69 


C 


1991 


1.00 


1.00 


1.00 


1994 


1.00 


1.00 1.00 


D 


1991 


1.67 


1.80 


1.71 









Table 7. The equivalent widths for the images in Q2237+0305, normalized with respect to measurements for image C at each epoch. 
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Figure 5. As for Figure g| presenting the spectra for image B and image C in 1991, which has been scaled so that the continua match. 
In the top panels, the B image is represented by a thick line, whereas the scaled C image is a thin line. The difference between these 
spectra is presented in the lower panel and the spectrum of the B image subtracted from the scaled spectrum of the C image. 
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Figure Al. Upper panel, model-fit (solid line) including con- 
tributions from the quasar and galaxy components (dashed lines) 
and the data (heavy solid line). X = 0.545. In the wings, 

at Pixels< 86 and Pixcls> 98 only the galaxy contributes to the 
model fit. The lower panel presents the residuals (data- model), 
in counts (solid line and left hand axis), and normalized by the 
noise (dashed line and right hand axis). 



Al Initial Reductions 

Bias subtraction and fiat-fielding were accomplished using 
standard procedures within the IRAF reduction package. To 
accomplish sky subtraction, the sky level was estimated di- 
rectly from the CCD frames by defining strips either side of 
the the quasar spectra. These strip were 15 pixels in extent, 
starting at least 30 pixels from the quasar spectra. At this 
distance the contribution of the extended galaxy is negligi- 
ble. The average of the median values of the sky in each strip 
was adopted as the sky value at each wavelength increment. 
There was no evidence of a gradient in the sky along the slit. 
For each object-frame a companion error-frame, recording 
the noise per pixel - including contributions from sky and 
object counts, CCD read-noise and fiat-fielding uncertain- 
ties - was retained. 

Wavelength and flux calibration proved straightfor- 
ward. Uncertainties in the wavelength calibration, derived 
from the residuals of individual arc lines about the cubic- 
spline fits employed were 0.2A in the blue and 0.05A in the 
red. Relative fluxes, to an accuracy of ~10% across the spec- 
tra, were obtained following repeat observations of standard 
stars. Given the narrow slits used and the systematic loss of 
light due to atmospheric dispersion it was decided that the 
final spectra should not be flux calibrated. 

Following sky subtraction, each CCD frame consisted of 
two quasar spectra superimposed upon the more extended 
lensing galaxy. Coefficients to perform wavelength and flux 
conversions were retained but not applied to the target 
frames, i.e. no rebinning or scaling that would compromise 
the independence of the noise from pixel to pixel was per- 
formed. The proximity of the quasar images to one-another, 
~ l"2, and the presence of the galaxy light precluded the 
use of conventional extraction procedures. Instead, a pro- 
cedure whereby the composite profile was modeled by three 



components, the two unresolved quasars and the galaxy, was 
employed. 

A2 Spectroscopic Extraction 

The key to the effective application of a multi-component 
model fitting procedure to estimate the contributions of the 
two quasar components and the galaxy at each wavelength 
increment is to reduce the number of free parameters to 
a minimum. In principle, the number of free parameters is 
large as three object profiles, their locations along the slit 
and their relative heights all need to be determined. How- 
ever, a significant reduction can be achieved by taking ad- 
vantage of the well-behaved and slowly varying properties 
of the spectra as a function of wavelength. 

Specifically, the centroids of the three components along 
the slit (pixel position on the Y-axis of the CCD frame) 
varies smoothly (see below) as a function of wavelength 
(pixel position on the X-axis of the CCD frame). An an- 
alytic fit to the location of the centroids of the components 
as a function of wavelength allows the three component cen- 
troid positions to be specified prior to the model fitting at 
each wavelength. 

The profiles of the two (unresolved) quasars components 
are identical, reducing the number of independent profiles to 
be determined from three to two. Furthermore, the variation 
in profile-shape with wavelength, arising from the combined 
effects of the wavelength-dependent atmospheric seeing and 
scattering in the spectrograph, is well-behaved, allowing the 
shape of the quasar profiles to be specified at each wave- 
length increment. 

Thus, implementation of the model-fitting procedure 
involves two stages: In the first, the systematic variation of 
the positions and shape of the component profiles is deter- 
mined and the shape of the component profiles estimated. In 
the second stage the contributions of the three components, 
essentially the relative intensity of the three component pro- 
files, at each wavelength (pixel) increment is calculated using 
a non-linear minimization routine. 

In the first stage of the model-fitting procedure, an un- 
weighted fit, neglecting any galaxy component, is used to 
establish the centroids of the point spread function of the 
quasar components. The location of the image centroids can 
then be represented by an analytic function, reducing the 
number of free parameters to three: the width of the seeing 
profile, and the normalization of each of the quasar compo- 
nents. The concluding stage of the procedure consists of a 
fully weighted model-fit, including a model for the galaxy 
contribution. Before considering the details of the quasar 
and galaxy profiles adopted, an overview of the entire pro- 
cedure is provided. The analysis proceeds as follows: 

(i) An unweighted fit of a model consisting of two iden- 
tical image profiles (the quasars) is undertaken for each 
wavelength increment in the CCD frame. Five free parame- 
ters, the profile centroids, the profile heights and the profile 
width (effective seeing) are determined via a Fortran90 im- 
plementation of a modified version of the AME0BA algorithm 
of Press et al. (1988). The initial location for the profile 
centroid (along the slit) is specified interactively for the first 
wavelength increment. For subsequent wavelengths, the lo- 
cation of the centroid from the previous wavelength incre- 
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Figure A2. Image centroids determined from the first stage of the fitting procedure and the analytic third-order polynomial represen- 
tations which allow the centroids to be specified at each wavelength increment. The data is from the red spectrum of images A+C, taken 
in 1994. 



merit is employed. The contribution from the galaxy profile 
is ignored at this stage. The unweighted fit is very sensitive 
to the peaks of the (quasar) point spread functions. 

(ii) An analytic function is then fit to the resulting pro- 
file centroid positions with the constraint that the profile 
separation remains constant (independent of wavelength). 
A third order p olyn omial representation produces excellent 
results (Figure A2). The zeroth-order terms were allowed 



to differ, while higher-order terms were forced to match for 
each of the quasars. 

(iii) A model fit, as in step (i) above, is then performed, 
but now variance-weighting is employed and the two image 
centroids derived from the polynomial representation deter- 
mined in step (ii) . Three free parameters are calculated, the 
profile width and the heights of the two quasar profiles. 

(iv) An analytic function is then fit to the resulting pro- 
file width as a function of wavelength. As for the centroid 
positions the dependence of the profile width on wavelength 
is smooth and well behaved allowing a very accurate repre- 
sentation using a simple analytic function. 

(v) Finally, a variance-weighted model fit is employed 
with the image centroids and profile width determined from 
their respective analytic representations. There are only 
three parameters remaining, the heights of the two quasar 
profiles and the galaxy profile. The relative contributions of 
each component from the fit averaged over one CCD expo- 



sure are illustrated in Figure Al 



An accurate representation of the quasar and galaxy pro- 
files and a quantitative measure of goodness-of-fit between 
data and model are required for the model-fitting approach 
to work effectively. 



p(r), consisting of a Gaussian core with exponential wings 
described the point spread function to very high accuracy. 
p(r) is defined: 



p{r) 



r < ri im 
r > rum 



(Al) 



where r is the distance from the centre of the pr ofile and ru m 
is the distance at which the two profiles join (Irwin 1985). 
By ensuring that the profile is continuous in the zeroth and 
first derivatives at ru m then 



p _ K2\o S N 

N 



V2 log N ' 



(A2) 



where is the fraction of the psf peak at which the profiles 
join. The condition that the profile is continuous in both the 
zeroth and first derivatives reduces the number of indepen- 
dent parameters to four; the central position, x c ; radial scale 
length, a; normalization, K; and the join fraction, N. 

Variations in atmospheric seeing produce changes in the 
PSF that are essentially confined to the central Gaussian 
core and the PSFs of all three standard star spectra could 
be modeled by fixing the join fraction, N equal to 35% of 
the peak profile intensity, thereby reducing the number of 
independent parameters necessary to describe the profile to 
three. 

As described in the previous section, the centroid and 
width (the a parameter) of the quasar profiles as a function 
of wavelength are determined prior to the final stage of the 
model fitting. Thus the profile height, or normalization, is 
the sole parameter to be determined for each of the two 
quasar spectra. 



A3 Model Components 

A3. 1 The Quasar Point Spread Functions 

The very high signal-to-noise ratio spectra of the standard 
stars, obtained using the identical spectrograph configura- 
tion on the same night as the observations of Q2237+0305, 
were used to define the form of the point spread function 
(PSF) for unresolved point sources. An analytic function, 



A3. 2 The Galaxy Profile 

The galaxy profile was modeled using a de Vaucoulers (1953) 
r« profile. The core radius was chosen to r e = 7"0, and the 
ellipticity to e = 0.31, in accordance with the observation of 
Racine (1991). The profile was smoothed with a PSF of the 
form de termin ed from the observations of the standard stars 
with the profile width defined by the value cal- 



fSection A3.1 



culated from the quasar spectra in the CCD frames. Virtual 
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4000 4500 5000 4000 4500 5000 
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Figure A3. The left panel shows a spectrum of NGC 3379, an 
E0 galaxy from the Kennicutt (1992) spectral atlas, redshifted to 
z = 0.039, the redshift of the lensing galaxy in the 2237+0305 
system. The right panel presents a galaxy spectrum extracted 
from the 1991 data for the B+D images. The spectrum has not 
been flux calibrated. A number of absorption features are common 
to both spectra, e.g.Call H&K AA3968, 3934, the 4000A break, 
the G-band A4300, H<5 A4102, H/3 A4861, and Mgl AA5173, 5175. 



slits of 0"7 and 0"91 were defined across the resulting surface 
brightness distribution at the same orientation as the actual 
observations. The light distribution was integrated across 
the slit at high resolution along the slit to give the spatial 
profile of galaxy light observed through the slit. The profile 
could be well reproduced using an analytical model consist- 
ing of two Gaussian components and a constant DC level, 
and this analytic representation was employed in the model- 
fitting procedure, with the height of the galaxy model profile 
as the free parameter. The normalization of the galaxy com- 
ponent and the two quasar components are obtained simul- 
taneously in the model-fitting procedure, thereby producing 
a galaxy spectrum in addition to the two quasar spectra. The 
galaxy spectrum extracted from one of the data frames, un- 
corrected for instrumental response, is presented in the right 
hand panel of Figure A3 . Fo r comparison, the spectrum of 
NGC 3379 ( |Kennicutt 1992| ), an E0 galaxy, redshifted to 
(z — 0.039) is shown in the left panel of the figure. The num- 



ber of absorption features in common between the spectra, 
including Call H&K AA3968, 3934, the 4000A break, the G- 
band A4300, U8 A4102, H/3 A4861, and Mgl AA5173, 5175, 
indicates the success of the fitting procedure in recovering 
a spectrum with the absorption line properties expected to 
be present in the galaxy core. 



A4 Goodness— of— Fit 

A standard X 2 (Aj) statistic was employed to ascertain the 
goodness-of-fit of the model to the data. The analysis proce- 
dure involves one-dimensional model fits, M(Xj,yi), to the 
observed distribution counts, S(Xj,yi), in the spatial direc- 
tion at fixed wavelength, (i.e. index Xj is fixed). 



[S(Xj,yi) — M(\j,yi)] 2 



Npi x ^ a' 2 (Xj,yi) 

i— 1 

where the model, M(Xj,yi) is: 
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Figure A5. The summed residuals from the model— fit to the 
synthetic data. Residuals (data-fit) in counts, in the spatial di- 
rection, averaged over the dispersion axis are shown as the solid 
line (left hand axis). The residuals, normalised by the noise at 
each pixel are shown as the dotted line (right hand axis). 



M{\j,Vi) = r + ^[Q 1 (X 3 ,y')+Q 2 (X ] ,y')+G(X J ,y')]d y ',(A4) 

where Qj and G are analytic models for the two quasar 
components and the galaxy respectively. These functions are 
integrated numerically over a pixel, where Ay is the pixel 
width. The noise, a 2 (Xj,yi), is: 



o 2 {X 3 , yi ) = S(X„ yi ) + Sky^+R 2 , 



(A5) 



where the first term represents the v N noise from the quasar 
and galaxy components, the second component the yHV 
noise from the sky, and the third component the read noise, 
R, from the CCD. 

The integrator used throughout was a Fortran90 im- 
plementation of the routine DQAG, a Gauss-Konrod general- 
purpose, globally adaptive integrator, available from the nu- 
merical mathematical software database GAMS This in- 
tegrator proved very capable in standard integral tests, es- 
pecially with the polynomial type integrands utilized in the 
evaluation of the x 2 statistic. 



A5 Fitting Procedure Applied to Synthetic Data 

To verify the effectiveness of the extraction technique, sim- 
ulated observations, closely mimicking the properties of the 
actual data, were generated and the ability of the fitting 
procedure to recover the (known) input parameters tested. 
The two quasar components, utilizing the composite quasar 
spectra derived by Francis (1991), were represented by PSFs 
with peak heights of 60 and 120 counts, and profile width 
a = l pixel. The components were separated by 3.7 pixels, 
with one centroid at pixel position 98.15, while the second 
was placed at 101.85. The galaxy component had a peak 
height of 15 counts and the spatial profile of Section A3. 2 



t http://gams.nist.gov/ 
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Figure A4. Results of applying the analysis to the synthetic quasar data. Lower panel shows the determinations of the quasar image 
centroids. Middle panel presents the quasar profile width measures. Top panel shows the quasar spectra recovered. The solid line in each 
panel represents the input values. 



Poisson noise from the object, a sky noise of 6 counts, com- 
bined with a read noise a n = 10 counts, were included. This 
resulted in a signal-to-noise ratio of ~2.5-3.5 for the sim- 
ulated quasar components, a value lower than the actual 
observations. 

Th e re sults of the fitting procedure are illustrated in 
Figure A4. The lower panel presents the measured image 
centroids from the first stage of the fitting procedure. Pi 
and P2, are the mean positions and the la scatter about 
the mean. The profile centroid separation is 3.69 ± 0.22 pix- 
els. The middle panel shows the recovery of the profile seeing 
width, a. The top panel presents the individual spectra re- 
covered by the analysis with the actual input spectra shown 
for comparison. Figure presents the corresponding mean 
residual (solid line and left-hand axis) in the spatial direc- 
tion. The dashed line represents the mean residual count 
normalized by the noise (right-hand axis) at each position. 
The residuals are small, typically ^ 15%, compared to the 
signal in the spectrum. The median value of the chi-squared 
merit function (Equation A3) over the fit is 1.20, a rea- 



sonable value considering the piece-wise and non-optimal 
nature of the fitting procedure. 

Recovery of the two quasar components from an artifi- 
cial data set with a signal-to-noise ratio and level of galaxy 
contamination closely matching that of the actual data gives 
confidence in the spectra obtained by applying the procedure 
to the observational data. 



A6 Robustness 

Differences will occur between the models employed in the 
spectral fitting, namely the point spread function and the 
surface brightness distribution of the galaxy, and the true 
functions. Such differences can introduce systematic errors 
into the extracted quasar spectra, affecting quantities such 
as the emission line equivalent widths. However, the parame- 
ters employed in the modeling are derived from observations 
of the lensing system and standard stars, and such system- 
atic errors should be small. These errors should, therefore, 
in no-way effect the overall conclusions of this paper. 
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